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High precision spectroscopy of muonic hydrogen has recently led to an anomaly in the Lamb shift, 
which has been parametrized in terms of a proton charge radius differing by seven standard devia¬ 
tions from the CODATA value. We show how this anomaly may be explained, within about a factor 
of three, in the framework of an effective Yukawian gravitational potential related to charged weak 
interactions, without additional free parameters with respect to the ones of the standard model. 

The residual discrepancy from the experimental result in this model should be attributable to the 
approximations introduced in the calculation, the uncertainty in the exact value of the Fermi scale 
relevant to the model and the lack of detailed knowledge on the gravitational radius of the proton. 

The latter cannot be inferred with electromagnetic probes due to the unknown gluonic contribu¬ 
tion to the proton mass distribution. In this context, we argue that muonic hydrogen acts like a 
microscopic gravimeter suitable for testing a possible scenario for the reciprocal morphing between 
macroscopic gravitation and weak interactions, with the latter seen as the quantum, microscopic 
counterpart of the former. 

PACS numbers: 04.50.Kd, 04.60.Be, 12.10.-g, 31.30.jr 


I. INTRODUCTION 


Precision studies of the hydrogen spectroscopy have played a major role in shaping our knowledge of the microscopic 
world and its understanding in terms of quantum mechanics and quantum field theory [I[ . Precision spectroscopy of 
hydrogen has now reached a point where the accuracy in the comparison to quantum electrodynamics is limited by the 
proton size, in the form of the root-mean square (rms) charge radius, Vp = To test quantum electrodynamics 

at the highest precision level, the proton size should then be determined with high precision from independent 
experiments. In the analysis of electron scattering experiments, a value of Vp = (0.895 ± 0.018) fm has been inferred 

0 : while higher precision determinations are possible using muonic hydrogen 0 . Since the more massive muon 
has a smaller Bohr radius and a more significant overlap with the proton, the correction due to the finite size of 
the latter is more significant than in usual hydrogen. However, a recent measurement 0 reported a value of Vp = 
(0.84184 ± 0.00067) fm, which differs by seven standard deviations from the CODATA 2010 value of (0.8775 ± 0.0051) 
fm, obtained by a combination of hydrogen spectroscopy and electron-proton scattering experiments. This anomaly 
corresponds to an excess of binding energy for the 2s state with respect to the 2p state equal to AE^^^p = 0.31 meV. 
Since we expect no difference between electron and muons in their electromagnetic behavior, due to the underlying 
assumed lepton universality for electromagnetic interactions, this has resulted in what is called “proton radius puzzle” 

ill. 

The existence of the proton radius puzzle seems confirmed by measurements of other energy levels allowing to 
determine the hyperfine structure with high precision 0, and has recently generated a variety of theoretical hypotheses, 
including some invoking new degrees of freedom beyond the standard model 0, [l0. In this context, some pioneering 
papers have already discussed high precision spectroscopy as a test of extra-dimensional physics for hydrogen PH, 
helium-like ions [10 . and muonium [l3 |. giving bounds on the number of extra-dimensions and their coup ling s. A 
recent attempt to explain the proton radius puzzle in the extra-dimensional setting has been discussed in [l5| . and 
values for the coupling constant necessary to fit the anomaly were inferred. 

The idea that extra-dimensions may be in principle tested with atomic physics tools is appealing also considering 
the paucity ofviable experimental scenarios to test quantum gravity [iM3- However, it would be most compelling to 
have a setting in which this may be achieved in the most economic fashion, that is, without necessarily introducing new 
free parameters conveniently chosen to accommodate a posteriori the experimental facts. In this paper, we provide 
such an approach by exploring the consequences of a tentative unification between gravitation and weak interactions 
already conjectured in (20|. In Section H, we introduce an effective potential energy between two pointlike masses 
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which recovers Newtonian gravity at large distances, while morphing into an inverse square law interaction with 
strength equal to the one of weak charged interactions at the Fermi scale, which coincides in our framework with 
the Planck scale. In Section III, we generalize this gravitational potential to the case of an extended structure like 
the proton, setting the stage for the evaluation of the Newtonian gravitational contribution to the Lamb shifts in 
perturbation theory. Section IV contains the main result of this paper, that is, the evaluation of the Lamb shift for 
the Yukawian component. The predicted contribution to the Lamb shift in muonic hydrogen is AE 2 s 2 p = 0.106 meV, 
to be compared with the experimentally determined value of = 0-31 meV, i.e. a factor 2.8 discrepancy. One 

potential source of discrepancy between our prediction and the experimental result is then discussed more in detail. 
This is then followed by the predictions for the expected contribution to the Lamb shift in muonic deuterium and 
a more qualitative discussion of the possible nature of the Yukawian potential of gravitoweak origin, including its 
selectivity towards the flavour of the fundamental fermions, this last feature being required in light of the manifest 
absence of a similar Lamb shift contribution for normal hydrogen. In Section V, we discuss possible tests in a 
purely leptonic system such as muonium, free from complications related to the extended structure of hadrons. We 
show that the contribution of the gravitoweak Yukawian potential is negligible with respect to the current precision 
achieved in the experimental determination of observables such as the Lamb shift and higher precision observables 
such as the absolute ls2s transition frequency. In the conclusions, we stress that a more accurate evaluation calls for 
the measurement of the gravitational radius of the proton, which is expected to significantly differ from the charge 
radius due to the gluonic energy density distribution for which no experimental access seems available. A qualitative 
discussion of other systems in which the effective Yukawian potential introduced here may give rise to observable 
effects, or may give significant constraints on its parameters, concludes the paper. 


II. AN EFFECTIVE POTENTIAL FOR GRAVITATION AT SHORT DISTANCES 

While we refer to [2^ for more details, we briefly recall here that the main idea we have explored is that what we call 
weak interactions, at least in their charged sector, should be considered as empirical manifestations of the quantized 
structure of gravity at or below the Fermi scale. This opens up a potential merging between weak interactions and 
gravity at the microscale, a possibility supported by earlier formal considerations on the physical consequences of 
the Einstein-Cartan theory |2l| . Various attempts have been made in the past to introduce gravitoweak unifications 
schemes , and a possible running of the Newtonian gravitational constant in purely four-dimensional models has 

been recently discussed [2l,|23. The conjecture discussed in [2^ relies upon identification of a quantitative relationship 
between the Fermi constant of weak interactions Gp and a renormalized Newtonian universal gravitational constant 
Gat, that is, we may write^ 


Gp = V2 Gn- ( 1 ) 

This expression holds provided that we choose Gn = 1.229 x IO^^Gat = 8.205 x 10^^ m^ kg“^ s“^. As an immediate 
benefit, the identification of the Fermi constant with a renormalized Newtonian universal constant via fundamental 
constants h and c allow to identify Fermi and Planck scales as identical, Ep = v, where Ep and v are, respectively, the 
renormalized Planck energy and the vacuum expectation value of the Higgs field, sometimes called the Fermi scale, 
avoiding then any hierarchy issue. As discussed in [^ . there are a number of possible tests of this conjecture that can 
span a wide range of energies, from the ones involved in the search for gravitational-like forces below the millimeter 
range [1^, to the ones explored at the Large Hadron Collider, with the spectroscopy of exotic atoms in between. 
Among the latter, we have outlined in the possibility that muonic hydrogen provides a suitable candidate, and 
here we make this proposition more concrete reporting an evaluation of this gravitational contribution to the Lamb 
shift in muonic hydrogen. 

A first tool required for our analysis is a proper interpolation between the two regimes of weak gravity at macro¬ 
scopic distances and the conjectured strong gravity/weak interactions at the microscale. The existence of Newtonian 
gravitation at large distances, with coupling strength given by the universal gravitational constant Gat, which can 
morph into weak interactions corresponding to a renormalized universal gravitational constant Gtv at small distances. 


^ Equation 1 differs from Eq. (2) in [Sfill since we have adopted in this paper a more rigorous definition of Planck mass as the one 
corresponding to the equality between the Compton wavelength and the Schwarzschild radius, i.e., h/{Mpc) = 2 GnMpIc^, the factor 
2 in the Schwarzschild radius having been omitted in the first analysis presented in 1201 . Equality of the Planck energy and the vacuum 
expectation value of the Higgs field, Ep = v, as in Eq. (4) of [2^1 requires the new prefactor in Eq. (1) of this paper. The value of the 
Planck length is invariant as for Ep = v we have Ap = hcjv, regardless of the definition of the Planck mass. 
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FIG. 1: Effective gravitational potential energy between an electron and a hypothetical pointlike proton versus their separation 
as from Eq. (2). The regions in which the inverse square law hold lie outside the two vertical dashed lines, as evidenced by the 
identical slopes. In these regions, we expect degeneracy for the 2s-2p contribution within perturbation theory, while in between 
the dependence on distance removes the degeneracy resulting in a mass-dependent contribution to the Lamb shift. Notice that 
this intermediate region extends over almost two decades starting at rmin — Ap and ending at rmax — 80Ap, therefore boosting 
by almost two orders of magnitude the influence of the lengthscale Ap at which quantum gravity effects are naturally expected 
to play a significant role. 


may be obtained by means of a generalized potential energy Ves for the gravitational interaction between two pointlike 
particles of mass mi and m 2 


14ff(r) = - 


GArmim2 




V Jr 

L / J 


( 2 ) 


wherein the intermediate expression we have introduced, as customary in the analysis of Yukawian components of 
gravity , the generic parameters a and A for the strength and range of the Yukawian component respectively, 

thereby specialized in the last expression to our case of interest, Q; = Gjv/GAr — 1 and A = Ap, with the renormalized 

Planck length Ap = 2hG]si/c^ = 8.014 x 10“^® m. 

Equation 2 reduces to ordinary gravity for r ^ Ap, whereas in the opposite regime of r <C Ap continues to have a 
1/r behaviour but with coupling strength proportional to Gfq. In both limits, the evaluation in perturbation theory 
of the average gravitational energy should give no difference since \/r potential is degenerate for states with the same 
principal quantum number and different angular momenta. Therefore, the difference we may evidence in this analysis 
will reflect the genuine deviation from a inverse square law characteristic of Yukawian potentials. As manifest in Fig. 
1 in the concrete example of an electron and a structureless, pointlike proton, the region in which the distance scaling 
of the corresponding force is not following the inverse square law extends over two decades starting at rmin ~ Ap. 


III. GRAVITATIONAL ENERGY BETWEEN A LEPTON AND AN EXTENDED PROTON: 

NEWTONIAN COMPONENT 

The evaluation of the gravitational energy between two particles is quantitatively different if one of them has an 
extended structure, as in the case of the proton. We schematize the proton as a spherical object with uniform mass 
ensity different from zero only within its electromagnetic radius Rp related to the rms charge radius - for a uniform 
charge density - through Rp = y'S/S r^, in such a way that Pp(r) = 2>mp/ {AttR^) for r < Rp and zero otherwise. 
With this assumption, the evaluation of the Newtonian potential energy - the first term in the righthand side of Eq. 
(2) - between the proton and a generic lepton of mass (£ = e, /x, r) yields 
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FIG. 2: Drawing (not to scale) of the relative proton and lepton coordinates for the evaluation of the Yukawian potential 
energy term due to the extended structure of the proton. The circle on the left side centered at point O represents the extent 
of the electric charge distribution of the proton, assumed uniform inside the sphere of radius Rp. The circle on the right side 
represents the range of the Yukawian interaction of a pointlike lepton located at point P. 
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( 3 ) 


The calculation of the energy contribution due to the Newtonian potential may be performed by means of standard 
time-independent perturbation theory applied to 2s and 2p states. There is no space anisotropy, so we will focus on 
the radial (normalized) components of the unperturbed wavefunctions which are, respectively 


R2sir) = 


(2a/!)3/2 y ae 


2 - 


e ; 


R2pir) = 


V3{2ae)^/‘^ ae 


r 

e . 


( 4 ) 


where ai is the Bohr radius, which also takes into account the reduced mass of the lepton-proton bound system. The 
calculation leads to 


{VNe)2s = {R2s\VNe\R2s) = - ——- X 


Rp 


16/32 


Fr 


-4K 


4^4 


/ + ^4 

0 / 16 


-4Fs 


-|- 4^2 




+ 00 


-4F2 


+ 00 


+ 4Fi 


+ 00 


(kjvf)2p = {R2p\VNe\R2p) = — 


G j^miTTip 

Rp 


48/3' 


;Fe 


P 1 

+ TT^4 
0 16 


P R +00 
+ —F 3 
0 24 ^ 


( 5 ) 


( 6 ) 


where we have introduced fi = Rp/at, and dxx'^e ^. 

a 

By straightforward algebraic manipulations, we obtain a more compact expression for the Newtonian potential 
energy difference between the 2s and the 2p state as 


^{VNe)2s2p = {VNe)2s — {VNe)2p 
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( 7 ) 
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FIG. 3: Plot of the Yukawian (continuous black line) and Newtonian (dashed red line) gravitational potential energies in 
hydrogen corresponding to Eqs. (3) and (10), respectively. The example assumes, to allow for an easier graphical comparison, 
interaction coupling strength and range as Gn = lO^Gjv, A = 2.0 x 10“^® m, and rms charge radius rp = 0.87 x 10“^® m, while 
the case considered in this paper has Gn = 1.229 x lO^^Gjv and A = Ap = 8.014 x 10“^® m, leading to a plateau value for the 
Yukawian component of the potential energy at small distance equal to ~ 0.52 /reV. The corresponding energy scale for muonic 
hydrogen is obtained by scaling up the vertical axis by the mass ratio of the two systems, leading to a plateau value ~ 0.108 
meV. 


from which it is easy to see that, in the (3^0 limit of a pointlike proton, A{VNe) 2 s 2 p —f GNmgmp(3^e~^/{2QRp) —>• 0. 
For a proton of hnite size, this contribution is different for hydrogen and muonic hydrogen since the gravitational 
mass mt appears both as a factor and in the expression for the Bohr radius at, upon which (3 depends. We notice also 
that this contribution tends to lift the energy of the 2s state, as this state is more sensitive to the finite size of the 
proton with respect to the 2p state, analogously to the case of the Coulombian attraction. However, it is also easy to 
check that the difference is absolutely negligible in our context, being about 37 orders of magnitude smaller than the 
experimentally observed anomaly in muonic hydrogen, so it cannot play a role in its understanding. This is consistent 
with a simple estimate of the gravitational contribution in muonic hydrogen based upon consideration of a pointlike 
proton: from the estimate for the absolute Newtonian potential energy Vg ~ GNm^rrip/aij, ~ 4.6 x 10“^"^ eV, and 
obviously the presence of an extended structure for the proton, and the differential energy taken in the Lamb shift 
evaluation, add up to further suppress the Newtonian contribution. This also indicates however that if the Newtonian 
gravitational constant is boosted by 33 orders of magnitude as in replacing Gn with Gn, the estimate for the absolute 
energy contribution is in the eV range, making its detailed evaluation relevant to the physics of muonic hydrogen. 


IV. GRAVITATIONAL ENERGY BETWEEN A LEPTON AND AN EXTENDED PROTON: 

YUKAWIAN COMPONENT 

The previous calculation sets the stage for the analysis of the Yukawian component, that is, the second term in 
the righthand side of Eq. (2). In a semiclassical picture, we expect that if the lepton is outside the proton by 
at least an amount ~ A, there will be a negligible Yukawian gravitational contribution. Likewise, once the lepton 
is completely inside the proton, there will be no Yukawian contribution either, as the uniform distribution of the 
proton mass will extert isotropic interactions. A nonzero value for the proton-lepton interaction occurs instead while 
the muon is partially penetrating inside the proton radius within a layer of order A. The exact calculation for the 
Yukawian potential felt by the lepton should proceed by integrating the Yukawian potential contributions due to 
each infinitesimal volume element inside the proton. Since this calculation involves exponential integral functions 
and it cannot be trivially performed analytically, we introduce two approximations. First, we truncate the Yukawian 
potential between two pointlike masses at distance r in such a way that 









6 


Vyiir) 


— cxG]\[ 

0 , 


miTTip 

r ’ 


0 < r < A, 
r > A. 


( 8 ) 


Therefore, the evaluation of the potential is carried out only in the region of intersection between two spheres, one 
of radius equal to Rp, the other of radius equal to the Yukawa range A. In using this approximation, the amplitude 
of the Yukawian potential is then overestimated in the 0 < r < A region, while it is underestimated in the region 
corresponding to r > A, as we will discuss more quantitatively at the end of the next section. Second, we further 
assume that the intersection region is a spherical cap, an approximation corresponding to consider Rp ^ A, which we 
will find a posteriori well satisfied. 

The potential felt by the lepton at a distance r from the proton center is then evaluated by integrating the 
infinitesimal potential energy contributions over the volume of the spherical cap of the sphere of radius A centered on 
the lepton location. With reference to Fig. 2, and using spherical coordinates for the infinitesimal volume dv, this 
can be written as 


/ dv ^27r n9 nX 

— = -aGNinipp / dcj) dcosO / ^d^, 

S J{) J J{r — Rp)/cos 9 


where cos0 = (r — Rp)/X. This leads to a Yukawian potential energy as follows 

0 < r < i?p — A, 
Rp — X < r < Rp, 
Rp ^ r Rp T A, 
r > Rp A, 


{ —2'KaGisi'nnpp}?, 

iraGNmf,pp[{r - RpY + 2A(r - Rp) - A^], 
-T:aG]simiPp[{r - Rp)'^ - 2A(r - Rp) + A^], 
0 , 


(9) 


( 10 ) 


in which continuity is ensured at all three boundary regions. As discussed above, this potential corresponds to a net 
attractive force only in the range Rp — X < r < Rp + A, and zero otherwise, and the plot of the corresponding potential 
energy is shown in Fig. 3. 

We evaluate the expectation value of the Yukawian potential energy in the 2s and 2p states, (VV^) 2 s and {VYe) 2 p 
according to hrst-order perturbation theory. The Yukawian component of the gravitational potential is still much 
smaller that the Coulombian potential even if it is coupled through Gn at short distances, since G^mtmp « 
e^/(47reo). The calculation leads to 
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( 11 ) 


( 12 ) 


where y = {Rp — X)/ai, z = {Rp + X)/ai, and v = X/ai. 

The difference between the two contributions finally yields, after lengthy algebraic simplifications, the expression 


HVyi)2s2p ^ {Vy)2s - {Vy)2p = "{^^(y) + [fl(y,/i) + KP,^)\] , (13) 
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FIG. 4: Dependence of the calculated gravitational energy difference between the 2s and 2p energy levels of muonic hydrogen 
AE 2 s 2 p versus the Yukawa range of short-distance gravity. We assume the effective coupling strength as Gat = 1.229 x IO^^Gat = 
8.205 X 10^^ m® kg“^ s“^ which gives rise to the Fermi coupling constant for charged weak interactions as discussed in [^ . The 
CODATA 2010 value of the proton rms charge radius is assumed for the gravitational proton radius. The experimental value 
for Aif 2 s 2 p = 0.31 meV is also shown (red cross), which is explained assuming a value of A = 1.35 x 10“^® m at the effective 
coupling strength and proton radius assumed above. The prediction of the model for A = Ap, corresponding to AE2s2p~0.106 
meV, is marked by the blue square. 


where f{y) = 2(2 — y)y^e ^ and 
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where G„ 


6 "e The sign of the correction is negative, i.e. AE 2 s 2 p < 0, indicating that the 2s state 


gets more bounded than the 2p state. This outcome is opposite to the case of the long-range Newtonian component 
in which the 2 s state gets a weaker binding since it explores more the proton interior where the net effect of gravity 
is smaller. This is due to the fact that the Yukawian potential does not fulfil the Gauss theorem and is marginally 
felt by the 2p state having small overlap with the proton region, unlike in the 2s state. We will indicate from now on 
the absolute value of AE 2 s 2 p, with the implicit understanding that it is negative. 

The result of this analysis is shown in Figs. 4 and 5, which summarize the main findings of this paper. In Fig. 4, 
we show AE 2 s 2 p as a function of the Yukawa range A for a fixed value of a corresponding to Gn as assumed in the 
text immediately after Eq. (1). For a value of A = Ap, we obtain a value of AE 2 s 2 p — 0.106 meV, about a factor 2.8 
smaller than the measured value. The experimental value of AE 2 s 2 p is instead obtained from Eq. (13) by assuming a 


value of A = 1.35 x 10 ^®m ~ 1.7Ap. In spite of the drastic assumption of uniform mass density for the proton inside 
its electric radius and of the two approximations used in the calculation, the proximity of the prediction on AE 2 s 2 p 
assuming first-principle parameters a = G^ jG^ — 1 and A = Ap as inspired by the conjecture in (^ . if not accidental, 
is rather remarkable. Considering the uncertainty in the relevant proton radius, we have complemented the analysis 



















FIG. 5: Robustness of the Yukawian gravitational potential with respect to the Yukawa range A and the rms radius Vp. The 
locus of the points in the (A,rp) plane which give rise to a 0.309 meV < AE 2 a 2 p < 0.311 meV is shown by the monotonically 
increasing curve. The range of values allowed, within one standard deviation, for the proton rms charge radius according to the 
CODATA value are delimited by the two vertical lines, the horizontal line corresponding to A = Ap. Complete agreement with 
the experimental value is obtained for instance by choosing A ~ 1.7Ap, while having Vp at its CODATA value, or by choosing 
the gravitational radius of the proton equal to ~ 0.7 times the CODATA value for the rms charge radius while keeping A = Ap. 


by showing, in Fig. 5, the combined effect of A and on the evaluation of AE 2 s 2 p- We have also evaluated the points 
in the (rp,X) plane which allow to obtain a value of AE 2 s 2 p in the (0.311,0.309) meV interval, see monotonically 
increasing line. The vertical band is the CODATA 2010 range of values for Vp, and the horizontal line corresponds to 
the value of A = Ap. To exactly verify our hypothesis we should have obtained a single intersection point, and this is 
not achieved within a relative accuracy conservatively estimated to be of order 80 % in A and 40 % in Vp. 

We further notice that the second approximation, consisting in considering the intersection region between the 
proton and the Yukawa sphere as a spherical cap, is well-satisfied since the proton radius is more than three orders 
of magnitude larger than the relevant Yukawa range. The truncation of the Yukawa potential expressed by Eq. ( 8 ) 
is instead a potential source of discrepancy, which can be removed by a numerical analysis with the full Yukawa 
potential. While this is certainly a subject for future investigation, it is worth to point out that the truncated 
potential over distance does, in general, overestimate the integral of the Yukawa potential. This can be seen more 
quantitatively for pointlike particles by considering the zeroth moment of the potential, its integral over the distance, 
Wy{t) = exp(—^/A)/^, and evaluating the coefficient a such that the zeroth moment of the truncated potential 

WTa{r) = equals Wy(r). Alternatively, one can evaluate a multiplicative coefficient 7 such that Wr-yir) = 

7 d^/^ equals lAV(r). In Fig. 6 we show the dependence of these two correction parameters on the distance r, 

evaluated in units of the Yukawa range A. It is evident that in the region of small r of interest for the evaluation of 
the integrals leading to the potential between the lepton and the proton as an extended structure both a and 7 are 
smaller than one. This means that in substituting the Yukawa potential with its truncated approximation one should 
either replace Ap —>■ crAp or AE 2 s 2 p —>■ lAE 2 s 2 p- Both procedures go in the direction of increasing the discrepancy 
discussed around Fig. 4, requiring a larger value of A or a smaller value of Rp to accommodate the experimental value 
of AE2s2p- 

All possible refinements of the calculations are limited by the fact that in this approach the mass density distribution 
is crucial, and we do not have any information on the density distribution for the most important component of the 
proton at the level of determining its mass, i.e. the gluons. This seems currently the most critical issue preventing a 
more quantitative comparison of the model with the experimental result. It is plausible that gluons, only sensitive to 
the attractive color interaction, tend to cluster more than valence quarks which are also sensitive to the electromagnetic 
interaction acting both attractively (between the quarks up and down) and repulsively (between the two quarks up), 
thereby decreasing the effective gravitational radius of the proton below its electromagnetic value. 

In spite of the limitation arising from the lack of knowledge of the gravitational proton radius, it is worth to proceed 
with the analysis of other systems in which this putative Yukawian potential may also give rise to observable effects 
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FIG. 6 : Correction factors a (black, continuous line) and 7 (red, dashed line) for the truncated potential in Eq. ( 8 ) versus 
the lower integration bound for the zeroth order momentum of the potential energy between two pointlike particles. Both 
corrections confirm that the overall contribution of the truncated potential is overestimating the Yukawian potential energy. 


and predictions. In this context, here we limit our attention to muonic deuterium and hydrogen. In particular, on the 
one hand we are able to make a firm prediction for the Lamb shift excess to be expected in muonic deuterium. On the 
other hand, the analysis of Lamb shift in hydrogen shows that there are potential troubles for the scenario envisaged 
here, which can be overcome only by relaxing the demand for an economic approach with minimal assumptions. 

Muonic deuterium [sol - ls^ may be another viable test of our hypothesis since we expect a significant difference due 
to the deuterium mass and radius, apart from complications due to the internal mass distribution of deuterium. The 
calculation for muonic deuterium proceeds by introducing the deuteron charge radius [ 3 ^ recommended by CODATA 
2010 as Td = (2.1424 ± 0.0021) fm, and obviously adjusting the mass as corresponding to the deuteron mass. We 
obtain, at A = Ap, an expected excess in the Lamb energy of 15.5 fieV, and of 42.6 /xeV at A = 1.35 x 10“^® m at 
which the Lamb shift for muonic hydrogen is fully explained. This corresponds to an anomalous frequency shift in the 
10 GHz range, smaller than the one in muonic hydrogen but still observable as well above the statistical (0.7 GHz) 
and systematic (0.3 GHz) uncertainties quoted in Q for muonic hydrogen spectroscopy. 

The predicted anomaly for the Lamb shift in (electronic) hydrogen is instead far too large, corresponding to 0.52 fieV 
at A = Ap and 1.47 fxeV at A = 1.35 x 10“^® m, if using the same GODATA 2010 proton rms charge radius, and this is 
potentially a big issue for the validation of the proposed model. The Lamb shift in hydrogen is known with a precision 
which cannot incorporate such a large contribution, corresponding to an anomalous frequency shift of 0.2 GHz against 
an absolute accuracy of 9 kHz, or 8.5 x 10“® relative accuracy. Among possible ways to circumvent this issue - next to 
consider the whole model as invalidated - is to assume that the effective interaction corresponding to the Yukawian part 
in Eq. (2) is flavour-dependent, and does not act among fundamental fermions belonging to the same generation. This 
is not dissimilar from recent attempts fl^ . to introduce new interactions which differentiates between leptons, 

and may be also related to the need to understand the role of the electron-muon universality in the theory-experiment 
discrepancy for the anomalous magnetic moment of the muon (STl - fs^ . A putative “off-diagonal” interaction naturally 
spoils the universality characteristic of gravitation. However, it should be considered that in a possible gravitoweak 
unification scheme the emerging structure should presumably incorporate features of both weak and gravitational 
interactions. The former is manifestly flavor-dependent, as shown in the presence of Cabibbo-Kobayashi-Maskawa 
(CKM) and Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrices for the charged current, so it is not a priori 
impossible that the interaction corresponding to the Yukawa component in Eq. (2) is highly selective in flavour 
content. This solution could inspire searches for models in which a “hidden sector” of the standard model includes 
intermediate bosons of mass in the range of the Higgs vacuum expectation value mediating interactions which have 
mixed features between the usual charged weak interactions and gravitation, for instance heavier relatives of the 
boson. As a final comment, we would like to also point out that no tests of the equivalence principle have ever been 
proposed involving other than the first generation of fundamental fermions. It would be interesting in this context to 
envisage tests of the equivalence principle, for instance free fall experiments involving strange matter such as neutral 
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K mesons. 


V. LAMB-SHIFT GRAVITATIONAL CONTRIBUTION IN PURELY LEPTONIC SYSTEMS 


We finally complement the analysis by evaluating the contribution of the putative potential in Eq. (2) for purely 
leptonic systems, such as muonium, a fi^e^ neutral bound state. Muonium is interesting from the theoretical viewpoint 
since, due to the pointlike structure of leptons, there are no complications arising from a mass distribution as in the 
proton case, while experimentally there are complications due to the finite lifetime of the muon. Unfortunately, as 
we will see below, the expected effect is way too small to be observed, but the related calculations shed also light on 
the reason for which muonic hydrogen is so effective at spotting a Yukawa contribution with range in the attometer 
scale. The calculation of the perturbation to the energy levels for muons due to a potential like the one in Eq. (2) 
leads to the following expressions 


(Ueff) 2 s = — 


(Vefl)2p = — 


GNmeTUfj, 

dOep 

Gnmenifj, 


4a 




1 + a 
1 + 


( 1 + 0 ^ 


( 1 + 0 ^ 


where ^ = Og^/A. The contribution to the Lamb shift is then 


4 ^ V 


(14) 

(15) 


A(14ff)2s2p = (f4ff)2s — (14ff)2p = —Ot — — -^- ■ 


2 

.( 1 +^ 


4,2' 

(TTeF^O+^J ’ 


(16) 


in which the Newtonian 1/r contribution is obviously canceled out since it is equal in the 2s and 2p states. In the 
limit of Oep A of interest in our considerations, i.e. ^ ^ 1, the leading order term is the first in the right-hand side 
so we get a simplified, yet accurate expression as 


(17) 

We obtain A(14ff)2s2p — 1-85 x 10“^^ eV, corresponding to a frequency shift of 4.48 x 10“® Hz. The precision in 
the determination of the Lamb shift in muonium, of order 1.5 % pol - ld^ . is not enough to observe the predicted 
effect in any foreseeable future. A precision observable in muonium is provided by the ls2s transition frequency. The 
contribution to the Is state due to the effective gravitational potential is 


(Ueff)ls = 


GNmetn^ 


1 -b 


Of 


( 1 +WJ 


The contribution to the ls2s frequency then is 


(18) 


( 19 ) 

2 ai^ 

which leads to a value of the same order of magnitude (just differing by a factor 7) of the 2s2p contribution evaluated 
above in Eq. (17). Although the ls2s transition frequency is experimentally determined with a precision of 9.8 
MHz (^ . the expected signal is still too low to provide an observable frequency shift. The situation in principle 
is less pessimistic for a state (the so-called “true muonium”), as the Lamb shift and the ls2s contributions 

get increased by a factor (m^/me)'^/8, leading to ~ 10 Hz for the 2s2p frequency shift, but the very experimental 
feasibility of muonium is still under debate [^, due to the complications of producing a bound state of two, rather 
than one, unstable particles. 

Nevertheless, this exercise shows the efficacy of muonic hydrogen with respect to purely leptonic, structureless, 
systems, in which the Bohr radius and the Yukawa range are the only relevant lengthscales. The presence of an 
intermediate lengthscale, the proton radius, gives a boost in the expected contribution of muonic hydrogen with respect 
to its muonium counterpart by an amount equal to {a^/Rp)^ ~ 8 x 10®, as arguable by means of simple dimensional 
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considerations. In other words, while leptonic systems get a contribution which is quite far from experimental 
observation, semileptonic systems such as muonic hydrogen instead get an amplified contribution due to the extended 
structure of the proton which shifts down the relevant lengthscale from ai to i?p, although at the price of introducing 
a parameter, the proton radius, which is not yet known from first principles. This presumably shows the existence 
of a narrow window for observing effects due to quantum gravity at the attometer (Fermi) scale. In any event, the 
fact that muonic hydrogen discriminates between 2s and 2p states, with respect to a Yukawian component of strength 
much larger than Newtonian gravity and acting on the attometer scale, should make it the first effective “quantum 
gravimeter”. 


VI. CONCLUSIONS 

In conclusion, we have shown that a solution to the “proton radius puzzle” is potentially available by means of an 
effective Yukawian potential originated by the morphing of Newtonian gravitation into weak interactions at the Fermi 
scale as conjectured in (^ . without additional parameters with respect to the ones already present in the standard 
model. Leaving aside the need for removing the approximations used in the evaluation of the Yukawa potential due 
to the proton as an extended object, either the relevant Yukawa range is not coinciding with the Planck scale, or the 
relevant proton radius is not the one determined using electric probes, or a combination of both factors. Among all 
these considerations, it is important to point out that the electric charge distribution of the proton investigated by 
using leptonic probes is not accurately representative of the mass distribution. Within the standard model gluons - 
which account for a large fraction of the proton mass - do not have electroweak interactions and therefore charged 
lepton or neutrino scattering off protons cannot provide any information on their spatial distribution. Furthermore, 
we have been forced to assume that the electron does not interact with the proton via the same effective Yukawian 
interaction as the muon, since the expected anomalous contribution in the former case is exceedingly large with respect 
to what is observed in the hydrogen Lamb shift. Although this spoils the requirement for simplicity and universality, 
it is plausible in light of the complex structure of weak interactions which are flavour- and generation-dependent. 

Lastly, it may be worth to briefly discuss possible consequences of our model for future experimental investigations. 
A very important proposal, as the one described in on muon-proton scattering, will not add insights on this issue, 
yet it will be critical in corroborating the existence of the anomaly and to pinpoint any possible deviation from the 
electron-muon universality. Unfortunately, the intrinsically short lifetime of the r-lepton of 291 fs makes it unfeasible 
to perform spectroscopic measurements on tauonic hydrogen, for which the predicted gravitational anomalous Lamb 
shift is even larger than that for muonic hydrogen. Also, provided that strong interaction effects can be effectively 
filtered out at the required accuracy level, a systematic study of precision spectroscopy of other exotic atoms with 
Bohr radii below the picometer range (see [1^ for a pioneering analysis), such as protonium [13,13, may allow to 
confirm or disproof our model in a reasonable time-frame. 
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